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Abstract 



We revisit our recently proposed formalism for dealing with the heavy quark 
limit of the quark Dyson- Schwinger equation. An ambiguity inherent in the 
original version of method is identified and resolved. Our reanalysis illustrates 
the importance of correctly accounting for the effect of hard gluon momenta 
on the heavy quark self energy in the vicinity of the bare fermion mass pole. 
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I. INTRODUCTION 



In a previous paper |3J] we have developed a formalism for dealing with the quark Dyson- 
Schwinger equation (DSE) in the heavy quark limit — > oo where m R is the renormalised 
heavy quark mass. The purpose of this brief report is to point out an ambiguity in the limit- 
ing procedure described in the previous paper, and to illustrate the procedure for extracting 
the correct heavy quark limit. The ambiguity arises from failing to take proper account of 
the contribution of hard gluon momenta to the loop integral in the DSE. We begin with a 
brief summary of the renormalised DSE. 



II. DYSON-SCHWINGER EQUATION 

Our starting point is the renormalised quark DSE |2|] 

e'(p,a) = z 1 ( / i 2 ,A 2 )^ J A ^ D , u (p-qh»s(q)r u ( q ,p) 



(2.1) 



where we have used a Euclidean metric in which timelike vectors satisfy p 2 = — ^Minkowski < 0> 
and for which {7^,7^} = 25^. We write the renormalised quark propagator in the form 



S(p,n) 



ry • pA(p 2 , /i 2 ) + B(p 2 , ji 2 ) 



i<y-pa v (p ,fi ) + cr s (p ,fi ) 



(2.2) 



where <Jv(p 2 ) — A(p 2 )/(p 2 A 2 + B 2 ) and <J S (p 2 ) = B(p 2 )/(p 2 A 2 + B 2 ). The unrenormalised 
self energy is written 



E'ip, K) = i 1 -p [A'(p 2 , A 2 ) - l] + B'(p 2 , A 2 ) 
where renormalised and bare quantities are related by 

A(p 2 ^ 2 ) = l + A'(p 2 ,A 2 )-A'^ 2 ,A 2 ), 



(2.3) 



(2.4) 



B{p 2 ^ 2 ) = m R ^ 2 ) + B'(p 2 , A 2 ) - B'(fi 2 ,A 2 ). (2.5) 
The renormalisation scale is set such that 

Sip)' 1 2 2 =t 1 -p + m R ( f i 2 ). (2.6) 

The set of equations ( |2.1| ) to ( |2.6| ) together with the 'abelian approximation' Z\ = Z 2 can be 
solved numerically for the propagator functions A and B once the renormalised quark-gluon 
vertex function T M , renormalised gluon propagator D^ u , renormalisation point (fi,mn(n)) 
and cutoff A are specified. 

In a general covariant gauge, the gluon propagator takes the form 

9 2 D»„(k) = L v - k -^\ A(k 2 ) + g 2 ^, (2.7) 



2 



where £ is the gauge fixing parameter. In this note we shall assume the bare vertex or 
"rainbow" approximation for the vertex function, namely ^(p, q) = 7^. 

Using Dirac trace identities to project out from Eqs. ( [2.1| ) to ( |2.3| ) a pair of coupled 
integral equations gives 



A'(p 2 ,A 2 ) = l + 



4Zi /- A d 4 q 



3p 2 J (2vr) 



p ■ q + 2 



p-{p-q)q-{p- q) 



p ■ q — 2 



p- (p - q)q • (p - g) 
(p - ?) 2 



(p- g) 5 



(p- # 



A 



(p - g) : 



(2.8) 



B>(p 2 ,A 2 ) = -Z 1 



(2^ 



3A 



(p-qf + £ 



(p - qf 



(2.9) 



III. THE HEAVY QUARK LIMIT 

The heavy quark limit is taken by assuming the expansions 

A(p 2 ^ 2 ) = 1 + E ^ ( f ' , B(pV) = m^/i 2 ) + Eb^.k), (3.1) 

where we have defined the new momentum variable K and renormalisaton point k by 

K = p2 + m «, K = » 2 + m «. (3.2) 
2im/j 2im R 

The change of independent variable p 2 — > K is illustrated in Fig. 1 of Ref. |IJ. Our aim is to 
solve for the quark propagator in the vicinity of p 2 = —m R , or equivalently, K = 0. From 
Eq. |3.1| we have 



m R a v (p 2 ) = a Q (K)+0[ — \, a s (p 2 ) — —<jq(K) + O ( J , (3.3) 

2 \m R J 2 \m R J 

where we have defined 

aQ{K) = U<Tm~y ^ J: {K) = V B {K)-'E A (K). (3.4) 



Note that Eq. (|J) is valid in any part of the complex momentum plane in which Eq. (|3.1|) 
holds, that is, there is no assumption that we are in the vicinity of K = 0. 

For the purposes of expanding the DSE in the inverse quark mass we introduce the 
momentum substitutions generic to heavy quark effective theory 

Pf, = imv^ + kfj,, q^ = imv^ + k'^ (3.5) 

where v^v^ = 1. Typically we take t> M = (0, 0, 0, 1), which entails K = k& + 0(1/ m R ) in the 
vicinity of p 2 = —m%. Note also that the result S(p,fi) = |(1 + j4)aQ(K, tt) + 0(l/m) for 
the dressed quark propagator is only a good approximation for \k^\ « m R . 
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In Ref. |J the heavy quark DSE is obtained by substituting Eqs.( |3.3| ) and ( |3.5|) into 
Eqs. (|2.8p and (|2.9|) , and then using Eqs. (|2.4j ) and fl2~5l) to obtain an integral equation for 
the heavy quark self energy T,(K). The procedure involves truncating a l/m# expansion 
of the integrand of the DSE. However, for heavy quarks, and for a model gluon propagator 
with a realistic ultraviolet tail, the DSE receives contributions from values of k'^ (defined 
in Eq. (|3.5|) ) which are not small compared with mjj. Because of this the 1/rriR expansion 
must be applied with some caution. 

To understand in detail what can go wrong with the procedure, consider the model gluon 
propagator used in ref. H, namely the smeared Frank and Roberts propagator 



A(k 2 ) = (2tt) 



[ m 2 d 
a 2 7r 2 



1 _ g-fc 2 /(4m?) 

k~ 2 



A IR (A; 2 ) + 



An 2 d 



(3.6) 



where d = 12/(33 — 2Nf), Nf = 3 is the number of light quark flavours, m t = 0.69 GeV is 
a parameter fitted to a range of calculated pion observables, and the Gaussian width a was 
chosen to be 0.5643 (GeV) 2 . As we have indicated, this propagator has the appearance of a 
Coulomb propagator beyond the scale of the parameter mj. Eq. (|2.8|) then becomes 



A'(p 2 ,A 2 ) 



4Zi f A d 4 q 



3p 2 
+4n 2 d 



(27T) 



p • q + 2 



p- - q)q -(p-q) 

(p-q) 2 



A 



IR 



(p-q) 2 a v (q 2 ) 



:i+OP- g+ 2 ( i-o p ^-y" g) 



w(q 2 ) 
(p - qf 



(3.7) 



where £' = g 2 ^/(47r 2 d). The procedure employed in Ref. of taking the heavy quark 
limit by substituting Eq. ( |3.5|) and keeping only leading order in l/m_R contributions to the 
integrand is certainly valid for the term containing Air. This is because the exponentially 
damped gluon propagator ensures that the factor in square brackets is only sampled for 
k'^ << m_R. For the remaining Coulomb part, however, one runs into problems. 
Consider for a moment the identity 







d 4 q 



p ■ q + 2 



p- (g - q)q -{p-q) 
(p-q) 2 



w{q 2 ) 
(p - qf 



(3.8) 



One easily checks that the integrand of Eq. ( |3.8|) is equal to 



\p\ kkv(g 2 ) d 

sin 2 9 dd 



sin 3 9 



(p - qf 



(3.9) 



where \p\ = (p ■ p) 1 ^ 2 , \q\ = (q • q) 1 ^ 2 and 9 is the angle between p^ and q^ in Euclidean 4- 
space. Clearly the angular part of the integral, namely / I(\p\ , \q\ ,9) sin 2 9d9, is identically 
zero. In Landau gauge, £' = 0, the second term in Eq. ( |3.7| ) should therefore integrate 
to zero because of cancellations of significant contributions from all directions in Euclidean 
momentum space. In terms of the Feynman diagram corresponding to the original Euclidean 
DSE, Eq. ( |2.1|) , this means that in the heavy quark limit hard momentum from the incoming 
quark line can flow through the gluon propagator line, significantly altering the direction 
of the internal quark line. Considerable cancellations occur in the integration over angles, 
resulting in a diminished net contribution to Eq. 
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Inserting the change of momentum variables Eq. ( |3.5|) into the last term in Eq. ( ft.7| ) we 
shall see below that these cancellations do not necessarily occur order by order in l/m R . In 
order to construct a viable formalism for extracting the heavy quark limit, it is therefore 
necessary to subtract off the angular oscillatory behaviour of the integrand before expanding 
in the inverse quark mass. 

To this end, we add to Eq. (|3.7|) an arbitrary multiple of the identity Eq. (|3.8|) to give 



A'(p 2 , A 2 ) = 1 + (infrared part) + 



167r 2 Zid 
3p 2 



dS H( P , q ) av{q2 



(2vr 



(p-q) 



4 : 



where 



H{p, q) = (1 + e + X)p ■ q{p - qf + 2(1 - £' + \)p ■ [p - q)q -(p-q), 



(3.10) 



(3.11) 



with A arbitrary. 

We can now make use of Eqs. (|3.3j ) and ( |3.5|) to extract from Eqs. (|2[ 
and ( |3 . 1 0| ) an integral equation for H(K,k). For example, in Landau gauge (£ - 
choosing = (0, 0, 0, K) (which is allowable to zeroth order in l/m R provided K 



D 

= 0), and 
<< m R ), 



we obtain 



E(K, k) 



rA d 4 k' 



(2tt) 4 ik' 4 + k' 2 /(2m R ) + «) 



|k' 


2 A IR 


\K - k' A ) 


2 + 


k' 


2 


[k - k' A y + 


k' 


2 



;i + a) 





k 


2 




(K-k' A y + 


k' 


2 





4n 2 d 



4) 2 + 



-(K 



(3.12) 



The A = case was solved numerically in Ref. JJJ. For arbitrary A, however, we find that 
the solutions are heavily dependent on the irrelevant parameter A, so something is clearly 
amiss. 

As pointed out above, and in the Appendix to Ref. |1], the integrand of this equation 
is sampled for values of k! which are not necessarily small compared with uir. The terms 
a s(q 2 )/(p — q) 2 an d cv(q 2 )/(p — q) 4 m the Coulomb parts of integrands of Eqs. ( |2.9| ) and 
( p.lO| ) cannot be the cause of the problem, since, as pointed out below Eq. ( [3. 3D , they are 
well represented by the heavy quark approximations ^crg/(/c — k') 2 and |ctq/[?b_r(/c — k') A ] 
throughout the complex momentum plane. The problem must therefore lie solely with trun- 
cation of the 1/rriR expansion of H(p, q) in the integrand of Eq. ( |3.10| ). By careful choice of 
the parameter A however, we find that we can subtract of the troublesome angular oscilla- 
tory behaviour hidden in the identity ( |3.8| ) to leave an integrand which can be accurately 
expanded in the inverse quark mass even for large values of k'. 

In the Appendix we show that, with k^ = (0, 0, 0, K), the parameter choice 



a = e/3 - 1 



implies 



H{p, q) 



k'l 



O 



(- 

\m R 



(3.13) 



(3.14) 



5 



We are interested in solving the DSE equation in the vicinity of the bare quark mass pole, 
that is, for \K\ « uir. In this case, as m^j — > oo, the non-leading terms in the I/thr 
expansion of H do not contribute to the integral in Eq. ( |3.10| ) and our procedure for obtaining 
the heavy quark limit is valid. For any other value of A however, the l/m,R expansion of H 
picks up 0{k' /uir) pieces which contribute to the ultraviolet part of the integral and cannot 
be neglected. This is the source of the discrepancy. 

In Fig. [I] is plotted the modulus of the heavy quark propagator obtained by solving 
Eq. (|3.12|) in Landau gauge with parameter value A = — 1, consistent with Eq. ( |3.13|) . For 
comparison, we also plot the solution from Ref. |TJ corresponding to A = 0. With the correct 
value of A, conjugate poles have moved considerably deeper into the complex momentum 
plane, and away from the imaginary K axis, or equivalently, away from the real, timelike p 2 
axis. Very little evidence is left of the original bare propagator mass pole at K = 0. 



IV. COMMENT ON EUCLIDEAN AND MINKOWSKI METRICS 

Our formalism presupposes the quark DSE to be formulated in Euclidean space. This is 
standard procedure in most nonperturbative treatments of confining field theory including 
DSE calculations 0. Opinions differ regarding the question of whether this procedure is 
simply a matter of practical expediency, or whether it reflects a deeper truth about the 
nature of confining theories. In its most extreme form [||, the position adopted is that the 
underlying quark degrees of freedom can only be meaningfully considered in Euclidean space. 
If we adopt a less extreme point of view, it may be instructive to explore whether the angular 
oscillatory function in the integrand of Eq. Q3.8p finds its analogue in a Minkowski formalism 
of the quark DSE. That is, does the quark self energy contain significant contributions from 
different directions in physical Minkowski space time which cancel? 

Consider the case where the Minkowski space vectors p^ and are both timelike. 
We define \p\ = (p^p^) 1 ^ 2 and \q\ = (q^q^) 1 ^ 2 - Choosing p^ = (\p\ , 0) and q^ = 
(r] \q\ cosh^, q \q\ sinh^) where rj = ±1 and 9 is real, the integrand of Eq. ( |3.8| ) becomes 

, n p-(p- q)q ■ (p - g) 

(p - q) 1 



o"y(g 2 ) rj \p\ \q\ <Jv(q 2 ) d ( sinh 3 ^ 



(p-q) 2 ~ sinh 2 # d6 \ \p\ 2 + \ q \ 2 - 2rj \p\ \q\ cosh 6 J ' 

(4.1) 

and the integration measure within the forward and backward light cones becomes 

d A q=Y< J o \q\ d\q\J o sinh 2 9d9 J dtt. (4.2) 

The analogue of the integrand over angles at fixed |g|, which is zero in the Euclidean for- 
malism, is now 

f°° ,„ d ( ri sinh 3 9 \ 

Y / d9— —„ ^ . 4.3 

^Jo d9 \\p\ 2 +\q\ 2 -2ri\p\\q\ cosh 9 J 

It is difficult to conceive of a regularisation procedure which could render such a severely 
divergent integral equal to zero. A similar analysis of the case where p^ is timelike and q^ 



6 



is spacelike also leads to a strongly divergent integral. We are led to the conclusion that 
the subtraction procedure described in the previous section is only relevant to the Euclidean 
formalism. 

V. CONCLUSIONS 

We have carried out a reanalysis of the quark DSE in the heavy mass limit m# — > oo. It 
is observed that, for a model gluon propagator with a realistic asymptotic ultraviolet tail, 
the procedure used in our previous analysis may not lead to a unique solution for the 
dressed heavy quark propagator in the vicinity of the bare fermion mass pole. The problem 
lies with failing to account adequately for the contribution to the DSE loop integral from 
the ultraviolet part of the gluon propagator. In essence, this means that in the Euclidean 
formalism, the momentum of an intermediate quark state cannot be assumed to be aligned 
closely with that of the external quark line in the heavy quark limit. 

We have solved this problem by adding to the integrand of the DSE a function which 
integrates to zero over angles in Euclidean space, but which cancels off higher order in I/tur 
contributions to the integrand that would otherwise contribute to the zeroth order quark 
self energy. We note that the analogous function in Minkowski space, when integrated over 
the two sheeted hyperboloid defined by a constant value for is a strongly divergent 

integral, suggesting that our procedure is unique to the Euclidean formalism. In fact, if 
such a procedure proved to be necessary and valid in Minkowski space, it would bring 
into question the notion that backward propagating intermediate states can be ignored in 
the heavy quark limit. In the absence of a viable nonperturbative treatment of quantum 
chromodynamics directly formulated in Minkowski space, it is difficult to place a physical 
interpretation on a procedure which appears to be simply a necessary mathematical step 
within the Euclidean formalism. 

The practical effect of the correction is to shift conjugate propagator poles present in 
our previous solution deeper into the complex momentum plane, and further from the real, 
timelike axis. Our solution for the heavy quark propagator can be regarded as the first 
step in determining the Bethe-Salpeter (BS) amplitude for Q-q mesons. In reference || it 
was found that, if a simple Gaussian gluon propagator which does not accurately model the 
ultraviolet gluon behaviour is used, the pole in the heavy quark propagator is only slightly 
shifted from its bare propagator position at K = 0. It was further shown that the presence of 
this pole prevented solutions to the ladder approximation Q-q BS equation. The removal of 
propagator poles to remote parts of the complex plane exhibited by our corrected solutions 
should overcome this difficulty. 

A calculation of the Isgur-Wise function using an impulse approximation involving DSE 
inspired light quark propagators and bare heavy quark propagators has been carried out in 
reference ||. In this calculation the bound state Q-q BS amplitude was modelled using a 
variety of Gaussian-like functional forms, each using a single fitted scale parameter. It was 
found that the calculated Isgur-Wise function is insensitive to which functional form used 
provided a suitable choice was made for the scale parameter. The appropriateness of the 
fitting procedure employed in reference || can be tested once the next step of modelling the 
BS equation with the corrected heavy quark propagator is completed. 
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APPENDIX: EXPANSION OF H 



Inserting Eq. ( |3.5|) into Eq. ( ft.llj ) we obtain 

H(p, q) = H 2 (p, q)m 2 R + H x (p, q)m R + H (p, q), 
where, after setting k^ = (0, 0, 0, K), 

H 2 = (1 + A) [-3{k' 4 - Kf - |kf] + e [{k' 4 - Kf - |k'| 21 
Ht = z(l + A) \3K 3 - 3K 2 k[ - Kk' 2 + 3k' 2 k[ - 2Kk' 4 



(A.l) 



J 2 



[-K 3 + K 2 k' 4 + 3Kk' 2 - k' 2 k 4 - 2Kk 
H = (1 + A) \3K 3 k' 4 + 3Kk' 2 k' 4 - 2K 2 k' 2 - AK 2 k' 4 \ + £ 



—K k 4 - Kk >2 k' A + 2K 2 k L 



The k l2 k' 4 term in Hi and the Kk' 2 k' 4 in H can be disposed of by choosing A = £/3 — 1, 
which reduces H (p, q) to 



H(p, q) 



4fm 



R n,/|2 



2iK K 2 



m R m Rj 



(A.2) 



For the region \K\ « ma, the final two terms can be ignored without affecting the solution 
of the DSE to zeroth order in I/uir. 



S 



REFERENCES 

[1] C. J. Burden, Phys. Rev. D 57, 276 (1998). 

[2] C. D. Roberts and A. G. Williams, Prog. Part, and Nucl. Phys. 33, 475 (1994). 

[3] C. J. Burden and D.-S. Liu, Phys. Rev. D 55, 367 (1997). 

[4] S. J. Stainsby, and. R. T. Cahill, Phys. Lett. A146 467 (1990). 

[5] M. A. Ivanov, et al, Phys. Lett. B 416, 29 (1998); Phys. Rev. C 57, 1991 (1998). 



9 



FIGURES 

FIG. 1. The modulus |ctq | of the heavy quark obtained by solving the heavy quark DSE 
Eq. (3.12) in rainbow approximation with a Landau gauge, gaussian smeared Frank and Roberts 
gluon propagator. The renormalised heavy quark mass is trr = 5 GeV, and remaining input 
parameter values are given in the text. The steeper surface with an obvious pole is the solution 
found in Ref. [|J corresponding to A = 0. The shallower surface is the solution corresponding to 
the correct parameter value, A = — 1. 
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